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Two roads are presently being followed in order to establish the existence of a liquid- 
gas phase transition in finite nuclear systems from nuclear reactions at high energy. The 
clean experiment of observing the thermodynamic properties of a finite number of nucleons 
in a container is presently only possible with the computer. Performed with advanced 
nuclear transport models, it has revealed the first-order character of the transition and 
allowed the extraction of the pertinent thermodynamic parameters. The validity of the 
applied theory is being confirmed by comparing its predictions for heavy-ion reactions 
with exclusive experiments. 

The second approach is experimentally more direct. Signals of the transition are searched 
for by analysing reaction data within the framework of thermodynamics of small systems. 
A variety of potential signals has been investigated and found to be qualitatively consistent 
with the expectations for the phase transition. Many of them are well reproduced with 
percolation models which places the nuclear fragmentation into the more general context 
of partitioning phenomena in finite systems. 

A wealth of new data on this subject has been obtained in recent experiments, some 
of them with a new generation of multi-detector devices aiming at higher resolutions, 
isotopic identification of the fragments, and the coincident detection of neutrons. Isotopic 
effects in multifragmentation were addressed quite intensively, with particular attention 
being given to their relation to the symmetry energy and its dependence on density. 



1. INTRODUCTION 

Multifragmentation reactions are often primarily seen and discussed in the context of 
the nuclear liquid-gas phase transition. Phase transitions in nuclear matter are indeed 
very special because they occur on temperature, pressure, and energy scales many orders 
of magnitude away from those familiar from ordinary matter. It is an interesting question 
in what form existing concepts for phase transitions can be applied to the phenomena gov- 
erned by the strong force. The observation or proof of the existence of a phase transition 
has been an important objective of heavy-ion physics for quite some time. 

There is little doubt that a liquid-gas phase transition should exist in extended nuclear 
matter |H2|3j . It is also known that the corresponding conditions of density and temper- 
ature can be reached in reactions between finite nuclei. The problems encountered there, 
however, are not only related to the small size of the reaction system, i.e. to the severely 
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limited number of constituents, but also to possible dynamical effects which may hinder 
the homogeneous population of the phase space. A third difficulty arises from the limits 
in our control of the experiments. The identification of the complex reaction processes is 
incomplete, and even the detected fraction of a multi-particle event is only known with 
finite precision. To sort data according to impact parameter, excitation energy or tem- 
perature is limited in accuracy by statistical and systematic errors. To work with the 
best possible coverage and resolution is thus even more essential, which has motivated a 
continuing process of improving the experimental possibilities and devices. 

More recently, a new direction for studying multifragmentation reactions has been de- 
rived from the importance of the symmetry term in the equation of state and of its density 
dependence for astrophysical applications. Supernova simulations or neutron star models 
require inputs for the nuclear equation of state at extreme values of density and asymme- 
try |4|5|6|7j . Isotopic effects in multifragmentation and other types of reactions have been 
shown to be sensitive to the symmetry energy coefficient, which may permit its study in 
the range of densities explored during the various stages of these collisions |8|9|imiT] . 

2. EXPERIMENTAL SITUATION 

Data obtained with 4-7r detection devices have been dominating the field for several 
years |12|13j . Besides solid-angle coverage, also the granularity and the dynamic range in 
particle type and energy are important parameters. A new standard for granularity has 
been set with the CHIMERA detector installed at the Laboratori Nazionali del Sud in 
Catania f^. It comprises 1192 individual Si-CsI(Tl) telescopes which are arranged in a 
ring structure providing highest granularity at forward angles [12]. 

The obvious need for the coincident detection of neutrons is known since long ago. In 
heavy reaction systems, a major part of the released excitation energy resides in the neu- 
tron channels, and calorimetry without neutron detection has to rely on assumptions. For 
isotopic studies, the neutron-to-proton ratio is a primary observable. Measurements of 
neutron multiplicities and energy spectra in coincidence with 4-7r charged-particle detec- 
tion have been performed with the NIMROD detector installed at Texas A&M University 
at College Station [15]. The data sets collected for selected reactions have permitted a 
comprehensive comparison with transport models. 

Apart from 4-7r devices, dedicated setups of high complexity have been developed for 
specific experiments. The measurement of neutron-neutron and proton-proton correlation 
functions has revealed isotopic effects in the space-time properties of the fragmenting 
source at breakup J7j. Fragment detection with mass identification is a prerequisite for 
studying isotopic effects in the fragment channels. For this purpose, new high-resolution 
arrays were constructed at Michigan State University jTH] and Indiana University |195. 

Possibilities for new experiments with existing detection systems have also been ex- 
plored. By transporting the INDRA multi-detector [20] to GSI, a new range of energies 
and reactions has become accessible for high resolution studies with A-ir coverage. Proton- 
rich secondary beams have been employed in experiments with the ALADIN spectrometer 
at GSI for the study of isotopic effects in the fragmentation of relativistic projectiles. First, 
preliminary, results from this latter experiment are now becoming available [21], some of 
the results obtained in the former campaign will be discussed further below. 
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Figure 1. Caloric-curve data from several experiments as compiled by Natowitz et al. 
|29j . Temperatures deduced from measured spectra of light charged particles and from 
double-isotope ratios of light fragments are shown as a function of the excitation energy 
deduced from the measured momentum transfer or from calorimetry. The Fermi-gas 
model predictions for parameters K = 8 and 13 are shown by the dotted and solid lines, 
respectively. The open circles represent the expectation for a hypothetical nucleon gas. 



3. THEORETICAL EXPERIMENTS 

The idea of studying the equilibrium dynamics of a finite number of nucleons in a 
container has first been realized by Schnack and Feldmeier several years ago |22]- The 
Fermionic Molecular Dynamics (FMD) model was used to propagate systems consisting 
of small numbers of nucleons over sufficiently long times, so as to allow them to establish 
equilibrium. With a weakly coupled thermometer, the temperature was measured as a 
function of the excitation energy of the system which was varied. A first-order transition 
from a Fermi liquid to a Fermi vapour within the confining oscillator potential is clearly 
observed. It appears as a natural consequence of the short-range repulsion and long-range 
attraction of the potential used for simulating the nuclear forces. 

Theoretical experiments of this type, based on the Antisymmetrized Molecular Dynam- 
ics (AMD) model, were conducted by Sugawa and Horiuchi and, very recently, by 
Furuta and Ono [21]. The outcome in these cases is very similar. The observed latent 
heat which is a function of the chosen volume or pressure indicates a first-order transition. 
The agreement reached by these studies for the location of the critical temperature near 
12 MeV in medium-heavy systems is remarkable. 

To validate these results, it is essential that the ability of the FMD and AMD models 
of realistically describing nuclear dynamics are thoroughly tested, and here considerable 
progress has been made. In particular, the AMD model is in the process of being in- 
tensively compared to reaction data, and very satisfactory agreement has been obtained 
|16|25j . Both, the AMD and the FMD, have also been developed into successful theories 
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Figure 2. Distributions in as a function of the transverse and longitudinal radius of the 
source for weighted (left panel) and projected (right panel) fragment-fragment correlation 
functions (5 < Z < 7) for ^^sxe + ''^^Sn (from Ref. gU). 



for ab-initio calculations of the structure of light nuclei |26|27j . The idea of deriving sup- 
port for the nuclear liquid-gas phase transition from studies of cold nuclei seems rather 
intriguing. 

4. THERMODYNAMICAL PARAMETERS 

Temperature measurements are useful for localizing the breakup process in the nuclear 
phase diagram [2H1- Here, a remarkable consistency has been reached recently. The sys- 
tematics compiled by Natowitz et al. [211 of experiments selected according to standard 
criteria, but allowing for different methods, is shown in Fig. ^ The combined set of tem- 
peratures exhibit an apparent Fermi-gas-like rise at excitation energies up to 3-4 MeV 
per nucleon, a very slow rate of temperature increase at higher energies, and some in- 
dication of a further rise at the very high excitation energies above 9 MeV per nucleon. 
Temperatures of around 6 MeV prevail in the range of excitation energies 5 to 8 MeV per 
nucleon at which multifragmentation is the dominant reaction channel for heavy systems. 
On the scale provided by the Molecular Dynamics studies, they place multifragmentation 
well below the critical point and into the coexistence region where homogeneous systems 
cannot exist in equilibrium. 

The apparent temperature spread (Fig. H} disappears to some part if the measurements 
are sorted according to the mass of the studied systems ^29j. The lower temperatures 
in the fragmentation domain (fa 6 MeV) are associated with heavy systems of mass 
A 200 while the higher temperatures (~ 8 MeV) are observed for the light systems 
with A < 100. The steady rise of temperature with excitation energy for the very small 
systems may even indicate that here the breakup occurs close to the critical point jSU] . 
This systematic dependence has been associated with the concept of the limiting tem- 
perature at which excited homogeneous systems become unbound [HII. According to 
these finite-temperature Hartree-Fock calculations, the limiting temperature should also 
strongly depend on the isotopic composition, a prediction to be tested experimentally. 
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The critical density of nuclear matter is approximately one third of the saturation den- 
sity [2]. For the statistical models of multifragmentation, the assumption of similarly low 
densities has, very early on, been found essential for correctly reproducing the observed 
fragment multiplicities |32|33I34] . Many experiments have shown that the kinetic-energy 
spectra of emitted fragments are inconsistent with emission from a normal-density com- 
pound nucleus |35l36l37|38j . More direct information on source properties in coordinate 
space is accessible with the technique of interferometry. Here, the shape of the source, 
the extension of statistical descriptions to non-spherical sources j^HI, and possible con- 
nections to transparency and mutual stopping in the approach phase of the collision jlHI 
have received particular attention recently. 

A new result obtained with a promising method of interferometry is shown in Fig. |21 
Two types of directional fragment-fragment correlation functions are used to determine 
the volume and the shape of the source at breakup following central collisions of ^^^Xe 
_l_ natgj^ gg ]\/[gY pgj^ nuclcou The distributions given in the figure represent 
the quality of their fitting, under the assumption of negligible emission times, with re- 
sults for model sources with different longitudinal and transverse extensions, Rz and Rxy 
Large radius parameters giving evidence for expansion are obtained with either technique. 
The correlation functions constructed from the projections of the relative velocities, in 
addition, indicate a longitudinally expanded shape with axis ratio ^1.7. 




Figure 3. Mean thermal excitation energy (full circles) and collective flow energy (open 
circles) at freeze-out, extracted by means of the MMMC-NS model, as a function of the 
incident energy Eq/A for central collisions of ^^^Xe + "^*Sn at 50 MeV per nucleon and 
^^''Au + ^^''Au at 60, 80 and 100 MeV per nucleon. The lines are linear fits to guide the 
eye (from Ref. |1S1)- 



A rapid expansion, following an initial buildup of pressure, will appear as a collective 
fragment motion. Collective radial flow refers to simultaneous transverse and longitudinal 
collective expansions which are observed in azimuthally inclusive experiments. Collec- 
tive flow in central collisions of heavy systems is observed above a threshold energy of 
f» 50 MeV per nucleon [Hlf-I^ij and starts to rise rapidly at higher energies. This range 
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Figure 4. Distribution of ^max versus 
Abound for ^^^Au on Au at 1000 MeV 
per nucleon [H^. Conventional fission 
events are removed . The sliadings fol- 
low a logarithmic scale. 




bond probability 

Figure 5. Bond percolation: relative 
magnitude of the largest cluster as a 
function of the bond probability for cu- 
bic lattices of 81 and 27000 sites. The 
critical bond probability in the infinite 
sytem is 0.249. 



of intermediate bombarding energies has been systematically covered in the set of experi- 
ments conducted with INDRA at GSI [l^l- The statistical analysis of these data yields an 
equilibrated excitation energy E* of the fragmenting source which rises rather slowly over 
the energy range up to 100 MeV per nucleon, and a rapidly rising collective-flow energy 
Ecow- At 100 MeV per nucleon incident energy, both components have about the same 
magnitude (Fig. E)). 

The successful description of these highly dynamical processes with statistical frag- 
mentation models, including a decoupled flow |43J, raises the question why the effects 
of the collective motion are not more clearly visible. Collective flow should affect the 
partitioning of the system and, in particular, the survival probability for heavier frag- 
ments |44|45|46I47|48I45] . The effect of the flow on the charge distributions, on the other 
hand, may be simulated by adapting the value of the thermalized energy in the model 
description |45l48j . It is, therefore, quite likely that the flow effect is implicitly included 
in the parameters of the statistical description. At moderate flow values, the changes are 
expected to be very small |47|48j . 

An alternative approach to the question of the coexistence of equilibrated partitions 
and collective motion has recently been presented by Campi et al. [50j. Using classical 
molecular dynamics calculations and specific clustering algorithms, these authors find 
fragments to be preformed at the beginning of the expansion stage when the temperature 
and density are still high. The fragment charge distributions, reflecting the equilibrium at 
this early stage when the flow is small, remain nearly unmodified down to the freeze-out 
density at which the flow has fully developed. 
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5. LARGEST FRAGMENT AS ORDER PARAMETER 

In the search for an experimentally accessible order parameter of the phase transition, as 
observed in reactions, the largest fragment of the partition has appeared as a promising 
choice. It may be identified with the part of the system in the liquid phase, and it is 
correlated with the mean density which is the natural order parameter of a liquid-gas 
phase transition. 

Statistical model calculations for nuclear multifragmentation show that the disappear- 
ance of the dominating fragment is associated with a maximum of the heat capacity 
which is the more strongly pronounced the larger the system [211 . The disappearance of 
the largest cluster, with the variation of a suitable control parameter, has been identified 
as a prominent signal also in fragmentations of other systems as, e.g., atomic hydrogen 
clusters and the largest cluster is an order parameter in percolation theory [33]. In 
finite percolation lattices, the disappearance of a dominant largest cluster proceeds rather 
smoothly and with obvious similarity to the nuclear experiment (Figs. I4I5|) . 

The success of percolation models in describing the observed partitioning of nuclear 
systems [55 56,,57j and the apparent critical behaviour ^H] do not necessarily identify the 
transition as of second order but rather show that first- and second-order phenomena may 
be compatible in small systems Studies of the fluctuation properties of the largest 
fragment [UO] or of the relative magnitude of the two or three largest fragments jUI] are thus 
not sufficient to establish the first-order character of the observed transition. Alternative 
methods have been presented [62, , ,63j which, however, require calorimetry on an event- 
by-event basis which cannot be performed without assumptions. Also methodological 
aspects are under debate |64|65j . Comparative studies of several of these signals are 
presently being performed by the INDRA collaboration |ti6|67j . 

6. NEW DIRECTIONS 

The experimental study of particle and fragment production with isotopic resolution 
has led to the identification of isoscaling, a phenomenon shown to be common to many 
different types of heavy ion reactions |8I9I68I6^ . It is observed by comparing product 
yields Yi from reactions which differ only in the isotopic composition of the projectiles 
or targets or both. Isoscaling refers to an exponential dependence of the measured yield 
ratios R21 {N, Z) on the neutron number and proton number Z of the detected products. 
The scaling expression 

R2i{N, Z) = Y^iN, Z)/Y^{N, Z) = C ■ exp{aN + pZ) (1) 

describes rather well the measured ratios over a wide range of complex particles and light 
fragments For illustration, the scaled isotopic ratios S{N) = R2i{N, Z)/exp{i3Z) for 
the reactions ^^C + ^^2,i24gj^ ggg ^^^^ qqq j^^y nucleon, studied with INDRA at 
GSI jTT] , are shown in Fig. The slope parameter a is found to decrease with increasing 
centrality of the reaction. 

Of particular interest is the connection of the isoscaling parameters with the symmetry- 
term -Esym = liA — 2Z)'^/A in the nuclear equation of state which has been consistently 
established with several methods |9|lU|7Uj . The coefficient 7 |72|73j is proportional to 
the isoscaling coefficient a according to aT ^ 47 • {Zf/Al — where T is the 
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Figure 6. Scaled isotopic ratios S{N) for 
12Q ^ ii2,i24g^ at E/A = 300 MeV and 

600 MeV for intervals of reduced impact 
parameter 6/&max, with "central" indi- 
cating h/h^s_y, < 0.4. The dashed lines 
are the results of exponential fits accord- 
ing to Eq. (0) (from Ref. [HJ). 



Figure 7. Isoscaling coefficient a (top), 
double-isotope temperatures TneLi (mid- 
dle) and resulting 7app (bottom) for the 
reactions + ii2,i24gj^ _ 

300 MeV (full symbols) and 600 MeV 
(open symbols), clS cl function of the cen- 
trality parameter 1 — 6/6max (from [7T|). 



temperature and Zi and Ai are the charges and mass numbers of the two systems at 
breakup. The results obtained for the ^^C + ^^^'^^^Sn reactions are summarized in Fig. [7| 
Using the above equation and assuming that the isotopic compositions are practically 
equal to those of the original targets, an apparent symmetry-term coefficient 7app was 
determined, i.e. without sequential decay corrections for a. The results are found to be 
close to the normal-density coefficient for peripheral collisions but drop to lower values at 
central impact parameters. If the corrections for sequential fragment decay after breakup 
are taken into account the resulting coefficient 7 for central collisions is even smaller [71j. 

This result as well as those obtained for other reactions |ll|74j may be considered as 
first steps within a program of studying the symmetry term far from saturation density. In 
the future, such experiments may aim at profiting from the existence of new radioactive- 
beam facilities. Projectiles with large neutron excess seem to be best suited to study 
dynamical effects of the symmetry force, and central collisions at high incident energy 
will have to be chosen in order to reach the larger than normal densities at which the 
symmetry energy is least well known |75|76j . 

Stimulating discussions with my colleagues at the GSI and within the ALADIN and 
INDRA-ALADIN Collaborations are gratefully acknowledged. In addition, I would like to 
thank R.T. de Souza, H. Feldmeier, T. Furuta, J.B. Natowitz, H. Oeschler, and A. Pagano 
for providing me with graphics for the talk and for helpful comments. 
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